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ABSTRACT 

We examine energy conversion from accreting pair plasma to outgoing Poynting flux 
by black hole rotation. Our approach is based on a two-fluid model consisting of 
collisionless pair plasma. The electric potential is not constant along magnetic field 
lines, unlike an ideal magnetohydrodynamics approximation. We show how and where 
longitudinal electric fields and toroidal magnetic fields are generated by the rotation, 
whereas they vanish everywhere for radial flow in a split monopole magnetic field 
in a Schwarzschild black hole. Outgoing electromagnetic power in a steady state is 
calculated by applying the WKB method to the perturbation equations for a small 
spin parameter. In our model, the luminosity has a peak in the vicinity of the black 
hole, but is damped toward the event horizon and infinity. The power at the peak 
is of the same order as that in the Blandford-Znajek process, although the physical 
mechanism is different. 
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1 INTRODUCTION 

The Blandford-Znajek (BZ) process is widely discussed as 
a promising mechanism for the powerful central engines in 
active galactic nuclei, microquasars, and gamma ray bursts. 
In their seminal paper, Blandford & Znajek (1977) showed 
outgoing energy flux from the event horizon on the assump¬ 
tion of a steady force-free magnetosphere around a slowly 
rotating Kerr black hole. Owing to the simplified situation, 
the electromagnetic extraction of the rotational energy could 
be analytically demonstrated. This remarkable process has 
been studied, partly because the underlying assumptions are 
doubted and partly because realistic astrophysical relevance 
is very important. 

Global steady-state force-free magnetospheres are con¬ 
structed by numerically solving the relativistic Grad- 
Shafranov equation, in which there are singular surfaces 
that make careful treatment necessary (e.g., Uzdensky 2004, 
2005; Contopoulos et al. 2013; Nathanail & Contopoulos 

2014) . Higher-order corrections to the original split 
monopole solution are calculated with respect to the 
black hole spin (e.g., Tanabe & Nagataki 2008; Pan & Yu 

2015) . The Grad-Shafranov equation for magnetohydro¬ 
dynamics (MHD) equilibrium is much more compli¬ 
cated and difficult to solve (e.g., Takahashi et al. 1990; 
Nitta et al. 1991; Beskin & Par’ev 1993). A time-dependent 
approach may thus be preferable, and general relativis- 
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tic MHD simulations provide very interesting models (e.g., 
Koide et al. 2002; van Putten & Levinson 2003; Komissarov 
2004b, 2005; McKinney 2006; Komissarov & Barkov 2009; 
McKinney et al. 2012). Recently, very complicated but more 
realistic configurations of the magnetic fields have been nu¬ 
merically treated. Of particular interest is the dynamical 
process in gamma ray bursts. The accretion of matter and 
jet emission in the vicinity of a central black hole can be ex¬ 
plored simultaneously. Sometimes, transient features are also 
exhibited in the numerical simulations. Outgoing flow from 
a system of a black hole coupled with surrounding magnetic 
fields has been tested, but it is rather difficult to determine 
the most important elements from numerical results. 

The mechanism for converting rotational energy to out¬ 
going electromagnetic flux resembles that of a pulsar. How¬ 
ever, the origin of electromotive force in a unipolar induction 
model is not established in the black hole magnetosphere. 
Possibilities include the event horizon (Thorne et al. (1986) 
and the critique by Punsly & Coroniti (1990)), ergosphere 
(Komissarov 2004a; Toma & Takahara 2014), and the pair 
creation surface (Beskin & Kuznetsova 2000; Okamoto 2012, 
2015) (see also, (Punsly 2008; Beskin 2010) and the refer¬ 
ences therein). When the ideal MHD condition holds in the 
entire magnetosphere, the angular velocity Dp of the mag¬ 
netic field is a function of magnetic flux, and this character¬ 
izes the electric potential difference between magnetic field 
lines. This value is crucial in the BZ mechanism, which works 
in only a certain range of flp. There remain the problems 
of where and how it is specified. Toma & Takahara (2014) 
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argued that a breakdown of the ideal MHD condition in the 
ergosphere is essential for giving rise to electromotive force. 
However, their argument is qualitative by an almost analytic 
treatment. It is impossible to discuss the origin of electro¬ 
motive force in the framework of the ideal MHD, and it is 
necessary to study it beyond the approximation level. 

In this paper, we do not assume ideal MHD conditions; 
we instead consider a two-component plasma consisting of 
positively and negatively charged particles, whose flows are 
governed by electromagnetic fields and gravity. Maxwell’s 
equations are solved with source terms of electric charge and 
current derived by the plasma motions, thus obtaining self- 
consistent solutions. Our approach allows us to explore the 
origin of the electromotive force and outgoing electromag¬ 
netic flux, if it exists. There are few works applying a two- 
fluid model to astrophysical situations, particularly to the 
formalism in black hole spacetime(Khanna 1998), and sta¬ 
tionary pulsar models(Kojima & Oogi 2009; Petrova 2015). 
There are large numbers involved in the model to connect 
microscopic to macroscopic sizes, and this fact is an obstacle 
to numerical calculation. 

This paper is organized as follows. Electromagnetic 
fields in a Kerr spacetime are discussed in Section 2 us¬ 
ing 3-1-1 formalism. We provide a brief review, because the 
equations in many papers assume the ideal MHD condition 
and those without the condition are needed here. We also 
discuss plasma flows interacting with electromagnetic helds 
in curved spacetime. The equations of stream functions are 
derived. In Section 3, we present a model to investigate how 
black hole spin modifies plasma flows and electromagnetic 
fields, resulting in outgoing energy flux. In a Schwarzschild 
spacetime, both the flow and magnetic field are radial. That 
is, the magnetic field has a split-monopole configuration and 
the electric field vanishes. To consider slow rotation of the 
back hole, we adopt the perturbation technique and explic¬ 
itly obtain the results. We present our conclusions in Sec¬ 
tion 4. 

Note that in the following we assume axial symmetry 
and stationarity in electromagnetic fields and plasma flows. 
We use units of c = G = 1. 


2 BASIC EQUATIONS 
2.1 Electromagnetic fields 

In this section, we briefly summarize the Maxwell equations 
in a black hole spacetime. The Kerr metric in the Boyer- 
Lindquist coordinate is given by 


ds^ = —a^dt^ + — uidt)^, (1) 


where 


2 P ^ 2 2,2 2 o A 2,2 oj/f 

a — , p = r + a cos o, /\ = r + a — 2Mr, 

S'* 

yi2 

zu = — sm fl, S =(r + a ) — a A sm y. 


2Mar 
E2 ■ 


( 2 ) 


We consider stationary and axially symmetric fields, 
and the electromagnetic vectors E and B refer to quantities 


measured by a locally non-rotating zero angular momen¬ 
tum observer (ZAMO). The four-velocity u'' of a ZAMO in 
Boyer-Lindquist coordinates is given by 


dx'^ 

dr 


-,--P 

a a 


, [/3’',/3®,/3'^] = [0,0, -ujm]. 


( 3 ) 


In this paper, the notation i denotes the component of a 
vector in an orthogonal basis. Using vector analysis in a 
3-dimensional curved space, the time-independent Maxwell 
equations are expressed as (e.g., Thorne et al. 1986) 

V-E = 4-Kps, (4) 


V • B = 0, 


( 5 ) 


V X (qE) = (/3 • V)B - {B ■ 


( 6 ) 


V X (qB) = 47raJ-(,5-V)E-b (E-V)/3. (7) 


It is convenient to introduce three scalar functions 
G{r,6), S{r,6), and ^{r,9) to express these fields. These re¬ 
spectively represent the magnetic flux, poloidal current flow, 
and electric potential^. The following forms satisfy eqs. (5) 
and (6): 


B = 


VG X 

zu 


-b 


s ^ 

azu 


(8) 


E = - - X B = -i(V'I> - cjVG). 

a a a 


Their components are explicitly written as 

[Bf, Bg, Bi] = 


1 S 

p OLZU 


( 9 ) 


( 10 ) 


\Ef, Eg, El] — 


A1/2 1 

-($,. -a;G,. ),-(cE>,e -LoG,g ), 0 

ap ap 


( 11 ) 


Since E^ — B^ = {{ujzu/a)^ — l)|VG/ci7|^, it is easily 
found that the electric fields dominate inside the ergoregion 
if the toroidal magnetic fields and longitudinal electric fields 
vanish (S' = $ = 0). Equivalently, the magnitude of ‘E x B’ 
drift velocity Vd exceeds the speed of light inside the er¬ 
goregion, judj = |E X Bj/B^ = jaitxj/aj > 1. This indicates a 
breakdown of the MHD condition (Toma & Takahara 2014). 
As a result, longitudinal electric fields arise to redistribute 
charge density. 

The poloidal component (r, 9) of eq. (7) is given by 


4:najp 


yg X 

ZU 


( 12 ) 


^ In the literature, symbols such as A, or / are used 

instead of G, and the electric potential is related by the ideal 
MHD condition. We do not assume this condition, so a new set 
of symbols is used here. 
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The poloidal current flows along a constant line S. The 
toroidal component (</>) of eq. (7) is 


T’GH-- - wG,r) + w,e(<l>,s - cuG,e)] = -Anamji, 

ap^ ^ 

(13) 


where 2? is a differential operator given by 


®C = -V.|^VG 


A1/2 




G, 


( O! 


VtuA^G 




Finally, Gauss’ law (4) is explicitly written as 
A1/2 r/Ai/2^, 


p^zu 


le 


-cuG.r) 


(14) 




Once the electromagnetic fields are known, the outgoing 
energy flux through a surface at a radius r is calculated as 
follows (see the Appendix): 


Pem(r) = - jt)d6d^ 

In the above we have presented a general form of sta¬ 
tionary and axially symmetric electromagnetic fields, which 
are described by three functions. If the ideal MHD con¬ 
dition E ■ B = Q holds, we have 4? = 4?(G), or V'F = 
f2F(G)VG, where flp represents the angular velocity of the 
magnetic held. Moreover, when the force-free approximation 
p^E -I- j X i? = 0 is used, the azimuthal component leads to 
S = S{G). The electromagnetic fields are described only by 
the magnetic function G (the Grad-Shafranov equation). 
The number of differential equations decreases, but there 
remains the problem of specifying their functional relations, 
!:!f(G) and 5' = 5(G). 


2.2 Outgoing energy flux from horizon 

To evaluate Pem in eq. (16) near a black hole horizon (at 
fh = M -I- V — a^), we need the behavior of the func¬ 
tions 4? and 5. For a while, we assume that the ideal MHD 
condition holds, that is, 'l>,e = Q.YG,e- Near the horizon, 
the function 5 is specified by imposing the so-called Zna- 
jek condition (Znajek 1978), = —Eg, which is written 

from eqs. (10),(11) as 

S=i^,e-u}HG,g)-, (17) 

P 

where luh = aj(rH) = a/{2MrE) is the angular velocity of 
the black hole. This condition is equivalent to the incoming 
electromagnetic wave (Thorne et al. 1986). By substituting 
these expressions into eq. (16), the outgoing electromagnetic 
power is 

Pem(rH) = —-Hf(Hf — t<JH) [ — {G,g)^d9. (18) 

2 Jrn P 

This shows that the power Pem(rH) is positive when 0 < 
Hf < t^H. This is the outgoing energy flux from a rotating 
black hole in the BZ process. A remaining problem is how the 
angular velocity Hf of the magnetic field lines or equivalently 
the longitudinal electric potential 4? is determined. 


2.3 Fluid 

We adopt a treatment in which the plasma is modeled as 
a two-component fluid. Each component, consisting of posi¬ 
tively or negatively charged particles, is described by a num¬ 
ber density n± and velocity w±, which denote the values 
measured by a ZAMO. The proper density n± measured in 
the fluid rest frame is related to n± and a Lorentz factor 
7 ± = (1 — (v±/c)^)~^G i 3 y — (^n/~f)±, where an abbre¬ 
viation (n/ 7 )± = n±l'^± is introduced. We assume that the 
positive particle has mass m and charge e, while the nega¬ 
tive one has mass m and charge —e. The charge density and 
electric current are given in terms of n± and v± as 

Pe = e(n+-n_), (19) 

j = e(n+w+— n_u_). (20) 

The continuity equation for each component in the axisym- 
metric stationary condition is 

0 = V • [an*'yvp)± = V • {a.nvji)±, (21) 

where the factor a comes from the relation between the de¬ 
terminant of 4-dimensional spacetime metrics and that of 
3-dimensional space This equation is satis¬ 

fied by introducing a stream function E± (r, 9) as 

a{nvp)± = —VP± X ei. (22) 

From this definition, the number density is solved as 

n± = (Qtu)“^(|VP|(i>?-I-u|)"^/^)± 

= (a^)-i(|VP|(l-7-"-i^|)-'/")±. (23) 

From eqs. (20) and (22), the current flow function 5 of 
eq. (12) can be given by 

5 = 47re(P+-P-). (24) 

The law of momentum conservation in a stationary ax¬ 
ially symmetric state (Khanna 1998) is 

^Vj(aTf J = pni±fli+P'ii5^j_±en±7±(P-l-F±x.B)i ±(Pcoi)i, 

(25) 

where the first term denotes gravitational acceleration with 
g = —Vina, the second is a gravito-magnetic term with 
Biij = the third is the electromagnetic force, and 

the last is the collision term. Here we consider the cold limit, 
so that thermal pressure is ignored, and the stress tensor 
is Tf = (mn* 7 ^u*ii-’)±, energy density pmi = (mn* 7 ^)±, 
and momentum flux 5m± = (mn* 7 ^u*)±. Moreover, we ne¬ 
glect the collision term. As Khanna (1998) discussed, near 
the black hole horizon, the electron collision time becomes 
longer than the dynamical timescale, so the collisionless ap¬ 
proximation may be valid under certain conditions. Using 
the collision rate (Spitzer 1962) and free-fall timescale 
tg, the product is Uctg ~ 10“®, where the electron number 
density rie ~ 10^^{M/10~^ cm“®, estimated 
from the accretion rate, and thermal velocity at T = 10^^ K 
are used. Collisions are ignored in the dilute approxima¬ 
tion. Thus, the equation of motion (25) for each component 
through the global electromagnetic and gravitational forces 
is reduced to 

[(Vx7F)xF+V7]± i = 'y±gi+Hij{'yv^)±±— \e + v±x B\ , 

771 L J i 
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(26) 

where the left-hand side is written in the vector form. From 
eq. (26), it is clear that there are two conserved quantities 
along each stream line, namely, the generalized angular mo¬ 
mentum J± and the Bernoulli integral K± , which are equiv¬ 
alent to Utf, and ut for each fluid component: 

J± = {-ajyv^)± ± ^G, (27) 

K± = (ujzu'yv^)± + ay± ± (28) 

These quantities depend on only the stream functions F±, 
and the spatial distributions are therefore determined by 
F±, which is specified at the injection point in our model. 
The component of eq. (26) perpendicular to the stream lines 
gives 

\an±zD^ J 

= aw‘^n± [k'± - -b mo) ^J±) ±^S, (29) 

where wo = ujwla, and J'± and K'± are derivatives of J± 
and K± with respect to F±. 

We have thus obtained a system of equations that gov¬ 
ern the electromagnetic field structure and plasma flows. 
Four partial differential equations (13), (15), (29) should be 
solved with two integrals (27), (28) and the number density 
(23) derived by the stream functions. These are reduced to 
those obtained for pulsar electrodynamics in a flat spacetime 
(Kojima & Oogi 2009) by setting M — a — 0. These equa¬ 
tions for G, $, and F± are interdependent in a nonlinear 
manner, so iterative methods are needed to self-consistently 
solve a set of these equations. For example, assume that 
functions G, th, and F± are known. The azimuthal velocity 
and Lorentz factor 7 are determined by the integrals in 
Eqs. (27) and (28). The number density is calculated from 
eq. (23). Thus, the source terms, namely, the toroidal cur¬ 
rent in eq. (13) and the charge density pe in eq. (15), 
are calculated from the fluid quantities of both species. The 
source terms and complicated coefficients in the equation of 
F± are also calculated. A new set of functions is solved for 
these source terms with appropriate boundary conditions. 
This procedure is repeated until convergence. 

3 MODEL DESCRIPTION 

3.1 Parameters and normalization 

We now describe a general framework for determining the 
electromagnetic fields and charged flows described in the 
previous section. Here, we discuss physical parameters in¬ 
volved in our system. 

There are two independent parameters. One is a di¬ 
mensionless gyrofrequency, x = ojbM = eBoM/m, where 
Bo is a typical magnetic field strength. Associated with Bo 
is a characteristic number density ric = Bo/( 47 reM). The 
number density is reduced to the so-called Goldreich-Julian 
density if the timescale 2M is replaced by the stellar angu¬ 
lar velocity The actual number density is normalized 

by Aric, where A represents the multiplicity of pair plasma. 



Figure 1. A schematic illustration for electromagnetic fields 
and plasma inflows around a rotating biack hole. The magnetic 
held is radial without toroidal component, and electric field van¬ 
ishes in the outer region. Black hole rotation affects the plasma 
flow, and a new electromagnetic structure is induced in the in¬ 
ner region. There is a current sheet on an equator to support the 
split-monople magnetic field. 


We can express other physically meaningful quantities us¬ 
ing these dimensionless parameters x A. The normal¬ 
ized plasma frequency k with number density Xric is given 
by = tOpM^ = 47 re^(Anc)M^/m = Ax- The ratio of A 
to X is written as fc = A/x = (1/4) x {2mXnc)/ {Bq/ 8n), 
and represents the ratio of the rest mass energy density 
of pairs to the electromagnetic energy density. These num¬ 
bers X and A are very large in astrophysical situations. 
A typical value of x is 10^^(i3o/kG)(M/10®MQ), relevant 
to active galactic nuclei powers oc (MBo)^. The amount 
of pair plasma is unclear, but we here estimate it from 
the accretion rate M. Using the electron number density 
rie ~ 1 O®(M/ 1 O“^M_b)(M/ 1 O®M 0 )“^ cm“® near the hori¬ 
zon, we have k? — 10^^{MMe){M/ lY"M q). The 
other parameters are calculated as A = k^/x ~ 10® and 
k = fx^ ~ 10~^. See also Beskin (2010) for estimates for 
microquasars and gamma ray bursts, for which A ~ 10^®- 
10^^. It is true that A,x,«^® 1, but these values should 

be regarded as an order estimate with large uncertainties. 
In particular, the ratio k = X/x may drastically change in 
the cases of matter-dominated or magnetically dominated 
flows. Indeed, the activation condition of the BZ mecha¬ 
nism in the MHD flows is approximately given by fc < 1, 
for which the Alfven speed exceeds the free-fall velocity at 
the ergosphere(Komissarov & Barkov 2009). Magnetization 
parameter corresponds to 1/k. 

We provide an explicit model of electromagnetic fields 
and plasma inflows in a hemisphere (0 < 0 < vr/2), depicted 
schematically in Fig. 1. We assume that the electromag¬ 
netic fields at rout M are described hy B = Bo{M/r)'^ef, 
B^ = 0 and E = 0, where Bo is a constant representing mag¬ 
netic field strength. This condition differs from the wind so¬ 
lution by Michel (1973) in which toroidal magnetic field and 
electric fields are given by B^ sin 0 /r. 

Our concern is how the parameter flp is determined, so 
that the condition Up = 0 is used at rout. Such a split 
monopole magnetic field may be formed by strong currents 
on an equatorial disk, by which upper and lower hemi¬ 
spheres are detached. The electromagnetic fields are ob¬ 
tained by the derivatives of G = G/{BoM^) = (1 — cosO), 
S = S/{BoM) = 0 and $ = ^/{BoM) = 0, where G, S and 
are normalized dimensionless quantities. 
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As for the pair plasma, we assume neutral flow falling 
along the radial magnetic fields at rout M, that is, —>■ 0 
and Q! 7 ± —>■ 1. The stream fnnctions F± of both compo¬ 
nents should coincide there, since F+ — F- = 5'/(47re) = 0 
in eq. (24). Like the magnetic fnnction G, both functions 
are chosen to be radial, F± — — cos0), where 

nc = {4:TTeM) is a characteristic number density and the 
minus sign denotes inflow vf < 0. We introduce dimension¬ 
less stream functions F± = F±/{XncM^) and the dimen¬ 
sionless number density h± = n±/(Anc). The relation (24) 
between S and F± becomes S = X{F+ — F-). 

Under these conditions at rout, the integrals J± and K± 
in eqs. (27) and (28) are explicitly given by 

J± = = TxMF±, K± = a7± = 1, (30) 

where the Lorentz factor at rout(i^ M) is for simplicity cho¬ 
sen as 7 ± = Q“^(« 1). Once J± is specified, the azimuthal 
velocity can be solved at any point from eq. (27) as 

vf =T^{F± + G). (31) 

'P ri77± 

The Lorentz factor 7 ± can also be solved from eq. (28) as 

a'y± = 1 ± X \iijM{F± -I- G) — $] . (32) 

Since we have K'^ = 0 and Jf = eq. (29) is 

reduced to 

(33) 

where uiq = utzojo.. 

In this model, we have imposed the conditions 5 = 0 
and $ = 0 at rout, so the electromagnetic Poynting power 
(16) is zero (Pem(rout) = 0). If the ideal MHD condition 
holds in rn < r < rout, then we have 4> = 0 everywhere, 
including the black hole horizon as the asymptotic limit, 
since the constant flp along any magnetic field line is zero. 
Consequently, no electromagnetic power is produced. Our 
concern is how power is produced in the presence of black 
hole spin. For this purpose, we have to consider non-ideal 
MHD effects. 


aw V ■ 


M^7± 

an±w'^ 


VP± = ±x 


3.2 Spherical flow 


Here, we discuss the structure of the electromagnetic field 
and plasma flows in Schwarzschild spacetime. An analytic 
solution is given in terms of dimensionless functions as 

G = -F± = 1 - cos 6», $ = 5 = 0. (34) 


The electromagnetic fields are explicitly written as 


> [Er, Eg, E^] = [0,0,0]. (35) 


The flow velocity, its Lorentz factor, and number density of 
the flow are given by 


[rf 1 I ] 


r ) 


7± = 70 = (1 - 
_ _ n± _ 1 

^ “ 71 .^r2(r -'2M)' 


(36) 


The number density n± appears to diverge at the horizon 
r —> 2M, but the proper density n*± = n±/ 7 o is finite every¬ 
where, as n±/(Anc7o) = M^/^/('\/2r®^^). The proper density 
and the magnetic field strength at r = 2M are n*± = Anc/4, 
and Bf = Bq/A = B^. Thus, Bo and Xric are reasonable val¬ 
ues near the black hole horizon, estimated in the previous 
section. 

This is a spherically symmetric solution, so that the 
conditions imposed at radius rout are retained everywhere. 
In particular, the electromagnetic power is zero everywhere. 
The energy flow Pq of matter across a sphere with radius r is 
obtained by twice the value in a hemisphere (0 < 0 < 7r/2): 


Po = —2 X A-KXricmM^ 


2X , .2 

- BoM^. 

X 


(37) 


Here, a minus sign denotes inflow, and the factor 2 comes 
from summing the contribution of positively and negatively 
charged fluids (see the Appendix). 


3.3 Effect of slow rotation 

We examine the effect of black hole spin on the spherical flow 
given by eqs. (34)-(36). It is rather difficult numerical work 
to obtain consistent solutions for G, $, and F± in eqs. (13), 
(15), and (33), as they consist of nonlinearly coupled par¬ 
tial differential equations. We here consider the rotation as 
a small parameter, and expand these functions as, for ex¬ 
ample, F± = XncM^{F -I- SF±). We limit our consideration 
to the first-order effect. The Poynting flux given by eq. (16) 
is a product of and 5. Both are zero for a = 0, but are 
modified by the first-order rotational effect. Poynting flux 
is thus produced within this approximation level. Within 
the first-order effect of the black hole spin, the only differ¬ 
ence from the Schwarzschild metrics is the function uj, which 
is approximated as cu = 2M^a*r“®, where a* = a/M is a 
dimensionless Kerr parameter. We assume a* > 0, which 
determines the direction of perturbed vectors. In this sec¬ 
tion, the symbols a and w denote those for a* = 0, that is, 
o? = \ — 2M/r and w = rsinS. 

We next develop the perturbation equations. We first 
note that the rotational effect in eq. (33) is the term wo = 
a~^ujw in only the lowest approximation, and that the term 
in the equation of the stream functions E± works in an oppo¬ 
site direction with respect to the fluid species. We therefore 
consider only a class of perturbations 5F+ = —SF-. From 
the perturbation of eq. (31), is expressed by 6F^ and 
SG, but after careful consideration we conclude that SG = 0 
and 6vf = Svr. We also have SG = 0 in the perturbation of 

!/> tj> 

eq. (13), consistent with 0 = 5j^ oc {Svf — Svr). From now 

on, we will use a function SF = SF+ = —SF-, and Svlr is 

</> 

given by 

Svh = —xM{w'yo)~^SF. (38) 

Under these conditions the perturbation of eq. (32) is re¬ 
duced to 

57 "*" = = TX7o'5‘h. (39) 

The acceleration, which is opposite that of the fluid 
species, originates from the perturbation of electric poten¬ 
tial. The perturbation of number density, eq. (23), is also 
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opposite to the species and is given by 


5n± 

no 


-^SFg+xho-^y^S^ 


(40) 


Finally, we consider the perturbation equations for 
54? and SF. General forms expanded with the Legendre 
polynomials Pi{0) are given by 54? = 'Y,hi{r)Pi(9) and 
SF = — + P)~^Pi{r)Pi,e{d)sm6. The slow rotation cor¬ 

responds to the dipole perturbation with Z = 1 , so that 
the components with I 7 ^ 1 are decoupled with the rota¬ 
tional perturbation. Therefore we have 54? = h{r) cos 9 and 
SF = (l/2)p(r) sin^ 0 Substituting their forms into the 
perturbation equations of Eqs. (15) and (33), after some 
manipulations we derive the following: 


d 

r 2dh\ 

r 2 2 

K, a 

2 ■ 

ds 

[ ds) 

. y2s 

s 2 


h + 


■\/2A 4a* 


(41) 


2 ^ ( dp\ 

ds \ ds) 





2s3/2 




‘2X0'* 

s5/2 • 


(42) 


Here, we use normalized radial coordinate s = r/M. We 
thus have a coupled set of second-order ordinary differen¬ 
tial equations. There are very large numbers involved in the 
hrst terms on the right-hand side, namely the squares of the 
plasma frequency k and gyrofrequency X- The last terms 
represent the effect of black hole spin a*. Our concern is the 
range x^ ^ s-nd = Ax 1, but k = A/x is not so large. 
Replacing = kx^ in eqs. (41) and (42), we neglect higher 
order terms of x~" {n > 2 ) except the derivative terms. 
Thus, eqs. (41) and (42) are approximated in forms suitable 
for WKB analysis: 



-2 ,, 




-X'^A 


s/i\ 
a ) 


]=X~"Js- 


(44) 


Here, s* in eqs. (43)-(44) denotes tortoise coordinate s* = 
r^jM— {r/M) + 21og(r/2M — 1) satisfying dst/ds = a~^. 
Other terms in eqs. (43)-(44) are 

U = {k/V2)a%-^/\ V = (^V2k - 

4 = J. = (45) 


The system of Eqs. (43)-(44) is rather simplified, since the 
source term in eq. (43), the higher order term ~ X~^i is 
neglected, and the potential terms JJ and V involve a pa¬ 
rameter k = A/x- 


3.4 WKB solutions 

We first consider the homogeneous solution of eqs. (43)- 
(44). We seek an approximate solution of the WKB form 
p oc exp(xlE(s)) and h oc exp(xlE(s)), where x~^(^ 1 ) 
is a small WKB parameter. Substituting them in, we find 
the leading-order solutions correct to order x~^ ■ The four 

^ The form SF allows stream lines to hit on the equator, and 
therefore poloidal current may go into and out it. This point will 
be discussed later. 


independent solutions (two pairs) given below are denoted 
by Two types are clearly distinguished among the 

three points: (1) oscillatory or growing/decaying behaviors, 
(2) the relative ordering between h and p and (3) their rel¬ 
ative sign. A pair of type I solutions is given by 

h'^ = eyip{±ix J U^^'^ds/) 

= exp(±ifi:is^^"^), (46) 

p± = -2x-'s'/= [2'/=A:(s + 2) - (47) 

where ki = (2^^'*/3) K and an overall constant from the inte¬ 
gral is adjusted in eq. (46). The solution represents spatial 
oscillation of the plasma, whose wavelength is ~ k~^M. The 
solution satishes the relations p^ ~ x~^ x in the 

large x limit, and h^p^ < 0 if the value in square brackets 
in eq. (47) is positive. The last condition is satished when 
k = A/x > 9-1 X 10~^. The functions and pf are regular 
toward the black hole horizon a —? 0, r —? 2M. 

Another pair of type H solutions is 

pj*; = exp(±x J V^^^ds/), (48) 

4 = 2^l\-^kcr^ [2^^\{s + 2) - pf,. (49) 

Here, we assume that V is positive everywhere. This con¬ 
dition is satisfied for the parameter k = \/x > kc, kc ~ 
2.3 X 10~^. When k < kc the potential V becomes negative 
in some range ri < r < r 2 , and the function in eq. (48) 
becomes oscillatory there. A whole solution is obtained by 
matching functions at ri and r 2 . We expect that such a so¬ 
lution is possible for only a discrete value of k, namely, an 
eigenvalue, and requires more careful treatment. In the fol¬ 
lowing, our consideration is limited to the case k > fee. Note 
that values in the square brackets in eqs. (47) and (49) are 
positive for that case as well. 

For later convenience, we approximate the integral in 
the exponential in eq. (48). With a constant K 2 = 2^^'^k, 
eq. (48) is reduced to 

pj*j « exp(±K.2S^^'*). (50) 

We numerically verihed that the approximation is good so 
long as K 1 and k > kc- 

We next discuss properties of the type H solutions 
in eqs. (48) and (49). These are exponentially grow¬ 
ing/decaying functions with a relation ~ x~^ ^Pu ^ Pn 
in the large x limit. This ordering is opposite to that in 
hf and p/^. As discussed in the next subsection, the sign of 
^nPn > 0 is critical for outgoing energy power. The func¬ 
tion oc a goes to zero at the horizon a —? 0, while in 
general diverges as hfj ~ a~^. The divergence in hu will be 
eliminated by appropriate combination of pj^, as discussed 
below. 

A general solution of eqs. (43) and (44) without the 
source term Js is expressed by a linear combination of 
four functions as h = 'Y,dth/[{s), and p = X)CnPn(s). 
The solution of the inhomogeneous equation is obtained 
by varying the coefficients c)/ as h = "Y, Cn{s)h//{s), and 
p = Y, ct {s)pt {s) ■ We put these forms into eqs. (43) and 
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(44), and find that dcf /ds^ oc x ^ s-nd that the functions 
satisfy the following equation of order x°: 

( 51 ) 

We neglect corrections to the type I solution, and consider 
only the effect of source term Js on cj*(. By integrating 
eq. (51), a particular solution of the inhomogeneous equation 
is written in a concise form as 

2 ^out 2 rs 

= ~ s3/4Vi/ 4 ^13/4^1/4 exp(-xl ds'\)di 

(52) 


/•out 

~ ~ 2i/4fci/2g3/8 £ a“'r""^®exp(-«2|si^"-Ci/"|)dC. 

(53) 

This is a method to solve inhomogenous equations in terms 
of a Green function constructed by the WKB approximation 
(see, e.g.. Bender & Orszag 1999). 

The solution p of eqs. (43) and (44) is in general a 
sum of pfi and pf, a growing solution pf is ignored since 
urout/M 3> 1. Even if small value p)') (rout) may be involved, 
the contribution exponentially decreases with the decrease 
of r. The regular solution both at the horizon and at infinity 
is given with a constant c)) by 

P^Pii + Cupjj. (54) 

The solution h (~ X~^p) is obtained by the relation (49), 
which may be unchanged within the lowest approxima¬ 
tion. Both pn and p)) go to zero toward the horizon, 
PiiiPn O’ whereas hu diverges in general as 

hu oc a~^. The regularity condition is 5Er = —~ 
and 5Eg = —(ar)“^((5<l?,6i — ojG^e) ^ in ZAMO vari¬ 
ables (Thorne et al. 1986). This means that the function hu 
should be finite at the horizon, and therefore the constant c)) 
may be uniquely determined as Cj") = —(pn/pO)a->o. That 
is, the divergence of /iii(oc a*’) is suppressed by an appro¬ 
priate of pii(oc a^). From now on, we consider this unique 
solution, which we call the type II solution for brevity. 

We next discuss the Znajek condition at the horizon. 
The condition of eq. (17) can be written in terms of the 
first-order perturbed functions as 

kp + X~^h = (55) 

where (54> = hcosO, 5S = 2X5F = Xpsiri^ 9, G — 1 — cos 9, 
SG = 0, and wh = a,/(AM) are used. This condition is a 
relation between p and X~^k, and can be easily incorporated 
in a type I WKB solution because their magnitudes are of 
order pf ~ x~^k.f. However, the ordering is opposite that of 
the type II solution, hf ~ x~^Pn- The leading-order WKB 
solutions are not satisfied with the condition (55), so that 
we have to take into account the higher-order WKB approx¬ 
imation as p oc exp(xIT) -s- p(i) -I- x“^P( 2 ) + X"^P( 3 ) H-, 

h —>• h(i) -I- X~^h( 2 ) + X~^h( 3 ) -!-•••, where the subscript (n) 
denotes approximation order n. The condition (55) provides 
a relation between p( 2 ) and h(i). In the type II solution, we 
approximate P(i) = 0 at the horizon, but have to treat a 
small correction X~^P( 2 )(<S 1) in a- more elaborate case. 



Figure 2. WKB solutions of p and h compared with numeri¬ 
cal integration of eqs. (41)-(42). The type II WKB solutions are 
plotted as solid curves, and the numerical results as dashed ones. 
The parameters are chosen as A = 1, x = 25 and a* = 1. 


3.5 Results 

We solve the differential eqs. (41) and (42) with the per¬ 
turbed Znajek condition of eq. (55) at the horizon and 
p = h = h^r = 0 at outer radius r/M = 10. The numer¬ 
ical solution is compared with the type II WKB solution 
p given by eq. (54) and the corresponding function h by 
eq. (49). Figure 2 shows good overall agreement except for 
the boundaries. The agreement in h is better. A minor dif¬ 
ference around r/M = 10 comes from the choice of outer 
boundary point; r/M = 10 is chosen in the numerical in¬ 
tegration, but it is infinity in the WKB solution. There is 
another difference only in p near the horizon. As discussed 
in the previous section, the function pn in leading-order 
WKB approximation goes to zero, but that of the numeri¬ 
cal integration goes to a finite value, satisfying the Znajek 
condition. Our WKB solution represents a piece of a numer¬ 
ical one, which involves another type I WKB solution and 
higher corrections of x~^- Note that the Znajek condition 
(55) leads to p « 0 at the horizon, if the function h is finite 
in large x limit. 

A set of differential eqs. (41) and (42) is not so compli¬ 
cated as to prevent integration. However, numerical results 
cannot be obtained with high precision for large dimension¬ 
less parameters because the characteristic lengths become 
very small. For example, we found that typical upper limits 
are x < 50 and A < 4. Big or small number is involved in 
numerical integration like e.g., exp(±Ks) ~ 10"*, 10“^ even 
for moderate case x ~ 10 and A ~ 1. Due to this limitation 
we use the type H leading-order WKB solution in place of 
numerical integration, and explore the behavior in a much 
larger x regime (x 1) relevant to astrophysical situations. 
This also provides an advantage of easily analyzing param¬ 
eter dependence. 

Figure 3 shows the electric potential d'h = hcos9{< 0) 
by the contours. The minimum is located on the polar axis 
at the horizon {9 = 0, r/M — >■ 2 or r^/M —>■ —oo), and 54? 
increases with either the 9 or the r* coordinate. The figure 
also shows that V5<1? is parallel to VG near the horizon; 
explicitly, 5<E>,0 = — (/i/M)G,e in the 9 direction. Therefore, 
the function h/M may be regarded as = —h/M, the 
angular velocity of the field line at the horizon. The velocity 
induced by first-order rotation is in the direction of black 
hole spin. A typical value of for the model shown in 
Fig. 2 is ~ O.OAat /M . By examining the behavior of h in 
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Figure 3. Contour of electric potential 0) in the 0-r* plane 

(left panel) and at the spherical coordinate (r, 6) (right panel). 
The minimum of (Sll? is located on the polar axis at the horizon 
(6 = 0, r = ch). 




0 


Figure 4. Contour of 5n+ /no{oc Spe) in the 6 -r* plane (left 
panel) and at the spherical coordinate (r, 9) (right panel). The 
minimum is located on the polar axis {9 = 0, r/M Ri 2.2). 

the WKB solution, we find that the value at the horizon is 
approximated as h « —a*/(2®= —a*/(2®'^^K). As the 
number density (i.e., A) increases, the deviation from the 
ideal MHD condition becomes small. Consequently, we have 
5$ —>■ 0 in the large k limit. 

We now consider the density perturbation induced by 
black hole rotation. Figure 4 shows contours of Sn+/no = 
—Sn-/no (oc 5pejno) in eq. (40). The minimum, negatively 
charged region for a* > 0, is located on the polar axis {6 = 
0). Unlike 54?, the minimum is not on the horizon, but at 
r/M « 2.2. The rate of charge density to the background 
number density goes to zero towards both radial directions 
r* —? ±00 and becomes neutral at the horizon and infinity. 

Figure 5 shows the stream function or current flow func¬ 
tion 25F = SS/X = psin^ 9 by contours. Since 5F < 0, both 
fluids are dragged in the rotational direction of the black 
hole, that is, Sv"/ > 0 (see eq. (38)). The effect is strong 
on the equator. The contours also show the poloidal current 
flow SS, where a half loop of current is formed around the 
minimum {6 — 'k/ 2, r/M « 2.5). In our treatment, we can 
never impose the condition SS^r = Xp^r = 0 on the equator. 


Figure 5. Contour of current function SS/X = 2dF{< 0) in the 
0-r* plane (left panel) and at the spherical coordinate (r, 9) (right 
panel). The minimum is located on the equator {9 = it/ 2,r/M 
2.5). 

This means that the poloidal current je oc —p,r is allowed 
to flow into or out from the disk. The function p^r changes 
sign at r/M « 2.5. The current emerges up [je < 0) from 
an outer point, say, ri/M > 2.5, flows along a constant line 
of 5S, and goes down to an inner point r\/M < 2.5 (see 
eq. (12)). However, there is no potential difference between 
ri and ri on the equator. Rather, the origin is attributed 
to the current sheet on the equator, which produces the ra¬ 
dial magnetic field in the vicinity of it. The right panel in 
(r, 0) coordinates shows a strong concentration of contour 
lines near the horizon. One might therefore assume a strong 
surface current on the horizon, but, as discussed in the pre¬ 
vious section, the function p in a type II solution goes to 
zero as ~ o? toward the horizon, so there is no poloidal 
surface current on the horizon. The left panel shows that 
the current vanishes for r* —? —oo. At the same time, the 
toroidal magnetic field SB^ ~ a also tends to vanish near 
the horizon. 

We calculate the electromagnetic energy power through 
a sphere of radius r. By integrating the perturbed expression 
of eq. (16) with respect to 6, we have 

5Pem(r) = ‘^XiBoMfhp. (56) 

The power 5Pem is positive for the type II solution, because 
hp > 0. In other words, outgoing flux is induced. In con¬ 
trast, the functions for the type I solution in which hp < 0 
represent incoming energy flux. 

Figure 6 shows the outgoing power SPem as a function of 
r for four models. The function increases with the decrease 
of r and has a maximum around r/M = 2.5, where the 
longitudinal electric fields and toroidal magnetic fields are 
significantly produced. However, the power declines toward 
the horizon and tends to zero owing to p —? 0. This energy 
flow is therefore produced by the black hole rotation, but 
is not related to the outgoing flux from the horizon. From 
numerical calculations, we found that 5Pem weakly depends 
on two parameters, (A, x) or ^ = ^/x)- 

Figure 6 shows that the peak slightly depends on k and 
increases with k~^, but does not depend on k. The math¬ 
ematical reason may be explained as follows: The WKB 
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10''3xPn 



Figure 6. Outgoing electromagnetic power SPem/ 

as a function of rjM. Four models are shown with k = 

0.025,0.04,0.1 and 1 from top to bottom. 

expressions eqs. (49) and (53) yield p ^ and h ~ 

so their product 5Pem ~ A x x 

~ in eq. (56), weakly depends on these parameters, de¬ 
spite p and h strongly depending on them in a different 
manner. The electromagnetic helds perturbations are ap¬ 
proximated as 5B^ ^ Ka,Bo and <5$ ~ K“^a*i3o, but the 
Poynting power, their product, is almost independent of the 
parameter k in the microscopic level. 

A typical value is of order 5Pem ~ 5 x [a^BoM)^ 
=8 X 10~^(a*PnA/)^, where the magnetic held strength 
Bn = Po/4 on the horizon is nsed. Compare this nnmber 
with the original estimate Pbz by Blandford-Znajek, calcn- 
lated with the ideal MHD and force-free approximations for 
a split monopole in the slow rotation limit. The result de¬ 
pends on the undetermined angular velocity flp of the mag¬ 
netic held lines. With the optimistic choice flp = <^h/ 2, we 
have Pbz = (a*PnM)^/6 ~ 1.6 x 10“^(a*PnM)^ in eq. (18). 
We therefore have SPem ~ 0.5Pbz. These are of roughly the 
same order due to the ambiguity of flp involved in Pbz- 
Finally, we calculate the conversion efficiency, which is a ra¬ 
tio of the electromagnetic power to that of the background 
inhows given by eq. (37): 

~ 2 X 10~^k~^al. (57) 

|Po| 

The efficiency increases with k~^, and reaches ~ 0.1a* for 
k ^ kc = 2.3 X 10“^. This mechanism is efficient in strongly 
magnetized how with less abundant matter, that is, smaller 
k. Contrarlly, the conversion becomes less active in high den¬ 
sity cases, and may not work for fc > 1. 


4 CONCLUSION 

We have explored the electromagnetic structure relevant to 
outgoing energy hux in a black hole magnetosphere. Our 
treatment was based on a two-huid description in which the 
ideal MHD condition is no longer assumed. By this formal¬ 
ism, we for the hrst time demonstrated how a longitudinal 
electric held is produced, even if it vanishes at the outer 
boundary. We hrst presented a general framework to con¬ 
struct a stationary axially symmetric structure of the elec¬ 
tromagnetic helds and pair plasma how. The formalism is 
an extension of our previous work (Kojima & Oogi 2009) 
for pulsar magnetospheres to a curved spacetime. 


For a dehnite result, we limited discussion to a slowly 
rotating Kerr black hole and used the perturbation method 
with respect to the spin parameter. We investigated how and 
where longitudinal electric helds and toroidal magnetic helds 
arise in the presence of black hole spin, although both helds 
are exactly zero in a spherical symmetric monopole model. 
To study the behavior at macroscopic lengths much larger 
than the plasma scales, we used a WKB method. By taking 
into account the hrst-order rotation, we found a unique solu¬ 
tion that describes zero helds at inhnity and hnite longitudi¬ 
nal electric potential at the horizon. This solution provides 
outgoing energy power, with magnitude on the same order 
as that of the original work of Blandford & Znajek (1977). 

In comparison with that work, we see a different physics 
involved in the plasma model. Split monopole magnetic 
helds as the lowest approximation and perturbations were 
used in both works. In both, the outgoing power originates 
from the black hole spin, but the origin is not the same. The 
black hole horizon, which is causally disconnected to the 
exteriors, plays a minor role; the energy power generated 
in our model is damped toward the horizon. The ergore- 
gion is often considered as the origin of longitudinal elec¬ 
tric helds and the outgoing power (e.g., Komissarov 2004a; 
Toma & Takahara 2014). The relation to the existence of 
the ergoregion is unclear within this work because we can 
not treat itiThe outer ergoregion radius coincides with the 
horizon, 2M in our lowest approximation of the spin. As 
discussed in Section 2, the ideal MHD condition should be 
inevitably violated inside the ergoregion for a purely poloidal 
magnetic held, but this argument may be too strong. The 
breakdown position is shifted outwardly in our pair model. 
Finally, Beskin & Kuznetsova (2000) discussed the impor¬ 
tance of the pair creation surface in determining the global 
MHD magnetosphere. The radius in a slowly rotating split- 
monopole magnetosphere is ~ 2.52M. Inflows and outhows 
are modeled as the boundary of both sides. In our model, 
the maxima of |(5pe|, \5B^\, and the luminosity peak are lo¬ 
cated around r ~ 2.2-2.5M. The pairs are assumed to exist 
outside the black hole, and to fall into the horizon in the 
background mean flows. We cannot treat outflows. It is im¬ 
portant to study the plasma behavior in the vicinity of the 
black hole horizon or ergosurface, which leads to outgoing 
power by the dragging. The crucial radius for the energy 
conversion is further pushed out by taking into account a 
realistic two-fluid Present result was derived under simpli- 
fled conditions such as steady state, split monopole, certain 
parameter ranges, boundary conditions and a leading-order 
WKB approximation. It is true that further study is needed 
to conflrm outgoing Poynting power generation by carefully 
examining each assumption. 
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APPENDIX A: ENERGY POWER 

In this appendix, we consider stationary energy flow. The 
energy momentum tensor is a sum of electromagnetic 
(Tern) ^nd matter parts. The outgoing energy per unit 

time through a constant radius r is obtained by integrating 
the energy conservation equation = 0. The 

power consists of a sum of the electromagnetic and matter 
parts: 

Pem(r) + P(+)(r) + P(_)(r) = - ^ {V^T[)d9d<j>. (Al) 

This total power is independent of r, but each part Pern, P(±) 
in general depends on r because of their interaction. The 
electromagnetic part Pem, that is, the outgoing Poynting 
flux, can be written as 

Pem(r) = - J (y^Te’n, 


where the energy momentum tensor is expressed by 
electromagnetic fields measured by a ZAMO, i = 
- tuwP")P‘^/47r. 

We next consider the plasma energy flow. Power in the 
+r direction is given for positively or negatively charged 
fluid components by 

P(±)(r) = - j,)ded^ 

= 27r {mK± =F e$)P±,sd6l, (A3) 

where ^ + uj'ujv^/ a)]± and Bernoulli inte¬ 

gral K± (28) are used. 

The total power eq. (Al), a sum of eqs. (A2) and (A3), 
is expressed as 

Pem + P(+)+P(-) f K±dF±, 

(A4) 

where we integrated by parts and used the relation S — 
4ne{F+ — P-). The first term on the right-hand side in 
eq. (A4) denotes the difference between the two boundary 
values with respect to 9, and is zero in appropriate cases, 
e.g., 5(0) = 4>(7r/2) = 0. The second term on the right-hand 
side is constant with respect to r, unless the plasma flows 
go out through boundaries 9a or 9),. We have thus confirmed 
that total energy flow is constant in a system consisting of 
an electromagnetic field and plasma. We consider the right- 
hand side in eq. (A4) for the linearized perturbations consid¬ 
ered in section 3. The first term is also zero up to the second 
order, since the boundary value of S^S + ^SS + S^SS be¬ 
comes zero at 0 = 0,7r/2. We assumed that K± = 1. and 
SF+ = —SF-, so that the sum in the second term becomes 
zero. 
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